Density functional theory for strongly interacting electrons 



0^ 
O 
O 
(N 



I 



I 

o 
o 



> 

p 

od 

o 
o 



Paola Gori-Giorgi/'^ Michael Seidl,'^ and G. Vignale* 
^ Laboratoire de Chimie Theorique, CNRS, Universite Pierre et Mane Curie, 4 Place Jussieu, 75252 Paris, France 
^ Afd. Theoretische Chemie, Vrije Universiteit, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands 
'^Institute of Theoretical Physics, University of Regenshurg, 93040 Regensburg, Germany 
'^Department of Physics and Astronomy, University of Missouri, Columbia, Missouri 65211, USA 

(Dated: September 15, 2009) 

We present an alternative to the Kohn-Sham formulation of density functional theory for the 
ground-state properties of strongly interacting electronic systems. The idea is to start from the 
limit of zero kinetic energy and systematically expand the universal energy functional of the density 
in powers of a "coupling constant" that controls the magnitude of the kinetic energy. The problem 
of minimizing the energy is reduced to the solution of a strictly correlated electron problem in the 
presence of an effective potential, which plays in our theory the same role that the Kohn-Sham 
potential plays in the traditional formulation. We discuss several schemes for approximating the 
energy functional, and report preliminary results for low-density quantum dots. 



Electronic systems are classified as weakly or strongly 
correlated depending on whether the potential energy is 
much smaller or much larger than the kinetic energy. In 
the first case the system is well described in terms of 
independent one-electron orbitals. In the second case it 
tends to crystallize in a "Wigner molecule" - a state in 
which the position of a single electron determines the 
position of all the others, even as any given electron is 
distributed in space according to the average density of 
the system. We will refer to this limit as the strictly 
correlated electron (SCE) limit. 

Which of the two descriptions is more accurate for a 
given physical system depends on the ratio of the aver- 
age inter-electron distance I to the effective Bohr radius 
a* = f? jra* e*"^ (including appropriately screened charge 
e* = e/e and effective mass m*): the strongly correlated 
regime occurs when a* 3> 1 and the weakly correlated 
regime when tja* <^ 1. For example, a system of N 
electrons trapped in the parabolic potential of a semi- 
conductor quantum dot ( "artificial atom" ) will enter the 
strongly correlated regime when the confinement length 
becomes larger than the effective Bohr radius. Similarly, 
a chain of hydrogen atoms becomes more strongly corre- 
lated as the distance between the protons increases. Fur- 
thermore all systems tend to become more strongly cor- 
related as the effective dimensionality is reduced. From 
a traditional quantum chemistry point of view, strongly 
correlated systems need very many (billions) of Slater de- 
terminants for a reasonable description of their physics. 
In other words, they are characterized by natural occu- 
pation numbers that are all close to zero. 

Strongly correlated systems pose a severe challenge to 
any many-body theory, because the electron-electron in- 
teraction cannot be treated perturbatively. In this paper 
we focus on the treatment of strongly correlated systems 
within the framework of density functional theory (DFT). 
DFT offers in principle a way to deal uniformly with both 
weakly and strongly interacting systems. In the original 
formulation of Hohenberg and Kohn (HK) the ground- 



state density and energy are obtained by minimizing with 
respect to the density p(r) the energy functional 



E[p]=F[p]+ / drv,^t{r)p{r) 



(1) 



where Dcxt(r) is the external potential and F[p] is a uni- 
versal functional of the density, defined as the expecta- 
tion value of the internal energy (kinetic energy T plus 
electron-electron interaction energy Vee) in the unique 
ground-state wave function that yields the density p(r). 

It is standard practice to carry out the minimization of 
Eq. (HD by the Kohn-Sham (KS) method ^^l. This method 
introduces the functional Ts [p] by minimizing the expec- 
tation value of T alone over all wavefunctions yielding the 
given p, and thus introducing a reference non-interacting 
system. The remaining parts of the exact energy func- 
tional are approximated. This works very well when the 
kinetic energy dominates, but runs into difficulties as the 
system becomes more strongly interacting. Indeed, a 
proper treatment of strong correlation is considered one 
of the two major problems facing DFT today (the other 
being a proper inclusion of van der Waals interactions). 

When the electron-electron repulsion dominates over 
the kinetic energy it could be much better to do the op- 
posite: define a model system in which one minimizes 
Vee alone over all wavefunctions yielding the given p, and 
approximate the remaining terms. This is precisely the 
approach we take in this paper. Thus, we present a rig- 
orous formulation of DFT that is alternative and com- 
plementary to the traditional KS-DFT approach in the 
sense that our reference system is a strictly-correlated 
system rather than the non-interacting KS system. This 
approach should be more suitable to treat systems whose 
density is not dominated by the quantum mechanical 
shells, but by incipient Wigner-crystallization effects. 

To formulate the problem precisely we start with the 
standard many-electron hamiltonian 



H = T+Ve, 



(2) 
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in which T = Si kinetic energy (we set 

2 

h = 1 throughout), Vee — z2i<j e\ri-r \ ^® ^^'^ clectron- 
electron interaction, and T4xt = X]i^cxt(ri) is the ex- 
ternal potential. The internal energy functional F[p\ is 
defined as [3] 



F[p]^mp]\f + V,,\^[p]) 



(3) 



where is the ground-state fermionic wave function 
uniquely associated with the density p. 

We further separate F[p] into a strictly correlated con- 
tribution T4e^^ [p] defined as the interaction energy func- 
tional for a system with zero kinetic energy (i.e., the 
minimum of Vee alone over all wavefunctions that yield 
the density p), and a remainder, which we call kinetic- 
decorrelation energy functional Ekd [p] '■ 



(4) 



Observe the close analogy between this separation and 
the conventional one involving the noninteracting kinetic 
energy, Hartree energy, and exchange-correlation energy. 
The functional Ekd[p\ corresponds to the expectation of 
the kinetic energy T plus the corrections to the expecta- 
tion of Vee- In this Letter, we explain how to construct 
the functional V"J,^^[p] and its functional derivative, we 
present three possible approximations for Ekd [p] , and we 
report preliminary results for low-density quantum dots. 

The functional V'ee^^[p] corresponds to the strong- 
interaction limit of the traditional adiabatic connection 
of DFT and was first addressed, in an approximate way, 
about 10 years ago 0, H, 0]. Only recently Seidl, Gori- 
Giorgi and Savin Q have found a general procedure for 
constructing the functional Vf^^ [p] , as well as other cor- 
responding observables such as the pair-correlation func- 
tion Q. Their construction can be viewed as the gen- 
eralization of the Wigner-correlated regime to any given 
non-uniform smooth density p(r). The key difference is 
that in the traditional Wigner regime (e.g., in Wigner 
atoms or molecules, or in the Wigner-crystal case), the 
electronic density is determined by the classical mini- 
mum of the hamiltonian without kinetic energy. Here, 
instead, we fix the density a priori, and for any given 
smooth (quantum mechanical) density we construct the 
corresponding V^^^ [p] . 

In the strong-interaction limit of DFT (no kinetic en- 
ergy, but fixed given density p) the admissible configura- 
tions of N electrons in d dimensions are restricted to a 
d-dimensional subspace of the full A'^d-dimensional con- 
figuration space 01 . We call this subspace A generic 
point of f7o has the form 



,(s) = (fi(s),....,fA.(s)) 



(5) 



where s is a rf-dimensional vector that determines the 
position of, say, electron "1", and fi(s) {i = 1,...,N, 
fi(s) = s) are the co-motion functions, which determine 



the position of the z-th electron in terms of s. The vari- 
able s itself is distributed according to the normalized 
density p(s) = p{s)/N. The co-motion functions are im- 
plicit functionals of the density, determined by a set of 
differential equations that ensure the invariance of the 
density under the coordinate transformation s — > fi(s), 
i.e. p(fi(s))c?fi(s) — p{s)ds 7|. They play the same role 
in our theory as the Kohn-Sham orbitals in the conven- 
tional formulation of DFT. In particular, the fj determine 
the functional VeJ~^^ [p] through the equation 



VeT^ip]^ I dsp{s) 



(6) 



just as the Kohn-Sham orbitals determine the non- 
interacting kinetic energy. Further, the total potential 
energy of a classical configuration 



£'pot(ri, ....vn) 



E 



e r,: - r,- 



Xl^scE[p](r»), (7) 



where usce [p] (r) is the external potential associated with 
the density p at zero kinetic energy, is constant on JIq 
and is expected to have a minimum with respect to vari- 
ations perpendicular to f^o, implying that its Hessian has 
d eigenvectors with null eigenvalue and N d — d eigenvec- 
tors with positive eigenvalue at every point on 0. 
The co-motion functions L have been constructed for a 
general spherical density [3, Q , while the solution in the 
general case is still the object of our on-going work. Al- 
though the functional T4e^^[p] depends on the density p 
in an implicit way through the co- motion functions (r) , 
its functional derivative with respect to p(r) is 



6p{r) 



(8) 



where the potential vsce [p] (r) satisfies the classical equi- 
librium equation Q 



V-yscE[p](r) = 



f.(r) 



N 

§k-f.(r)|3- 



(9) 



Combining Eqs. ([4]) and ([T]) we see that the minimization 
of the energy leads to the variational condition 



Sp{v) 



t(i-)+M- 



SEkdjp] 
Spiv) 



(10) 



where p is a constant (chemical potential). Assuming 
that we know how to approximate Ekd[p] and how to 
calculate its functional derivative, 



Wfed[p](r) 



5Ekd[p] 
Sp{r) 



we therefore reduce the energy minimization to the prob- 
lem of solving a strictly correlated system in an effec- 
tive potential wscE[p](r) = Woxt(r) + Vkd[p]{r)- An ad- 
vantage with respect to the KS approach is that the 
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co-motion functions fi(r) can be directly constructed 
from the density by integrating the differential equation 
p{ii{s))dti{s) ~ p{s)ds Thus, the simplest way to 
solve the SCE equations is probably by directly minimiz- 
ing the energy with respect to the density, using a proper 
basis set or a grid. 

The central problem is to obtain an explicit expres- 
sion for Ekd[p]- Here we discuss three approximations: 
a "first-order" approximation, a proper generalization of 
the local-density approximation, and the combination of 
the two. To define a "first-order" approximation, we par- 
allel the standard adiabatic connection approach to the 
calculation of the exchange-correlation energy. Namely, 
we introduce a fictitious hamiltonian 

Ho,=af + V,e + Vc. (12) 

where a is a positive constant that takes values between 
and 1 and Va — ^Q(i"i) is an external potential chosen 
in such a way as to yield the desired ground-state den- 
sity p(r) for every value of a. Notice that, at variance 
with the standard approach, we make an adiabatic con- 
nection between the reference system (strictly correlated 
electrons at a = 0) and the physical system (a — 1) by 
gradually turning on T rather than T4e- Making use of 
the Hellman-Feynman theorem it is easy to show that 

Ekdip] = I da{^^[p]\f\^^[p\), (13) 
Jo 

where ^a[p\ is the ground-state wave function associ- 
ated with the density p at coupling constant a. We 
have recently shown [9] that the exact expansion of 
Tq[p] = (*Q[p]|T|5'a[p]) for a ^ can be obtained from 
a classical harmonic analysis, which yields 

T„[p] =a-'/'rzp[p] + 0(a°), (14) 

where T'zp[p] is the kinetic energy associated with zero- 
point oscillations about the SCE solution. Inserting this 
expansion into Eq. yields a "first-order" , or zero- 
point (ZP) expression for Ekd[p]^ 

/Nd-d , 
dspis) Yl (15) 
n=l 

where u;„(s) are the Nd—d zero-point frequencies around 
the degenerate SCE minimum [9|. An example of cal- 
culation of Tzp[p] for spherical atoms from He to Ne 
is reported in Ref. 0. While Eq. is formally valid 
only in the limit a —> it is essential to appreciate that, 
in a physical sense, a small value of a is one for which 
a ^ £/a*. It follows that for a strongly correlated sys- 
tem, in which i/a* ^ 1, the physical value of a = 1 is 
already in the strongly interacting regime. 

A simpler approximation is a generalization of LDA, 

ElF^ip]- J drp{r)eMr)), (16) 
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TABLE L Relative % errors on the total energy of a model 
2D quantum dot consisting of two electrons confined in an 
harmonic potential ticxt (r) = ^uj^r'^. Columns as follows: KS- 
LDA are the results for standard Kohn-Sham LDA, SCE are 
the results obtained by setting Ekd[p] = in Eq. Q, SCE-ZP 
are the results obtained from Eq. p5|l . SCE-LDA are those 
obtained by using E^^^lp] of Eqs. and SCE-ZP- 

LDA are those obtained with Eqs. (fTg)l - ([Tg)) . 

where ekdip) is the kinetic-decorrelation energy of a uni- 
form electron gas with density p, given by 

(^kd{rs) ^ ts{rs) + exc{rs) - —■ (17) 

Here is the density parameter, defined by r^a* — 
{pBdy^/'^, where Bd is the "volume" of the d- 
dimensional ball of unit radius [lo| . The coefficient au 
determines the Madelung energy. Notice, again, the anal- 
ogy with the standard KS-LDA: in the Kohn-Sham for- 
malism, tg is treated exactly via the functional Ts[p], so 
that it is subtracted from the total energy of the elec- 
tron gas. In our case, Vgg^^[p] is treated exactly via 
the construction of the co-motion functions, and so the 
corresponding value for the electron gas (the Madelung 
energy) is subtracted from its total energy. Another way 
to look at Eq. (fT7|) is that LDA is uniquely defined as 
the approximation that makes Eq. Q exact for a uni- 
form density. As in standard KS theory, LDA has the 
advantage of being an explicit functional of the density, 
so that its functional derivative is easily calculated. 

A third approximation can be obtained by combining 
the zero-point energy and the LDA, 

^fc^r^'^'W = 2rzp[p] + j drp(r)6fc,_zp(p(r)), (18) 
where 

e-kd-zp{rs) = ts{rs) + e-xc{rs) - — - (19) 

is the total energy of the uniform electron gas from which 
we have now subtracted also the zero-point low-density 
term. Again, the spirit is to define LDA as the approxi- 
mation that makes Eq. (|18p exact for a uniform density. 

As a preliminary test for our construction and approx- 
imations we have used a simple two-dimensional (2D) 
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quantum-dot model consisting of two electrons confined 
in an harmonic potential Uext(r) = ^w^r^, with w mea- 
sured in effective Hartree*. At this first stag e we have 



used as inputs the exact densities of Ref. [11|. The 2D 



LDA functional is from Ref. [1^. The % errors on the total 
energy are reported in Table U where we compare, for dif- 
ferent values of w, the standard Kohn-Sham LDA results 
(column KS-LDA) with those from the SCE construction 
with Ekd[p] — (SCE), those from the "first-order" ap- 
proximation (SCE-ZP), those obtained by using E^^-^[p] 
(SCE-LDA), and those from the combination of the two 
(SCE-ZP-LDA). We see that when the system is weakly 
correlated (large uj) the KS-LDA result is superior, as 
the physics of the system is well captured by the non- 
interacting reference system. But as uj is lowered and 
the system becomes more correlated, the SCE construc- 
tion with its approximations for Ekd[p] becomes much 
more accurate than KS-LDA. The "first-order" (SCE- 
ZP) approximation gives the best results in the strongly- 
correlated regime (small oj), but the simpler SCE-LDA 
result is also accurate over a broad range of uj values [2l| , 
reducing the KS-LDA errors of a factor '--^ 4 up to '--^ 10 as 
the system becomes more correlated. The combination of 
"first-order" and LDA (SCE-ZP-LDA), instead, is always 
worse than the simpler SCE-LDA approximation. This is 
very similar to what happens in standard KS-DFT when 
we combine the exact first-order approximation (which, 
in this case, is the exact exchange) with the LDA cor- 
relation energy: the results are worse than when using 
LDA for both exchange and correlation. The poor per- 
formance of KS-LDA for small uj is due to the fact that 
the single Slater determinant is a very bad approxima- 
tion in this regime: exact-exchange yields much worse 
results, e.g., overestimating the total energy by ^ 60% 
at w = 2.211 X 10-3. 

We note that the two functional KI^^Ip] ^-nd Tzp[p] 
are formally independent of the statistics of the particles: 
that is because in the strongly correlated limit the par- 
ticles are distinguished by their relative positions, and 
exchange effects are in a first approximation negligible, 

— 1/4 

entering, formally, at orders e~" (i.e., the order of 
magnitude of the overlap of different gaussians in the 
zero-point oscillations). However, such exchange effects 
are in principle contained in Ekd[p]'- its LDA approxima- 
tion, for example, takes into account, in an approximate 
way, the fermionic nature of the system. 

In conclusion, we have presented a new formalism with 
high potential to treat strongly interacting systems such 
as low-density nanodevices. A representative of the class 
of problems that can be tackled by this formalism is 
the calculation of the addition energy of quantum dots 
[isl Aside from the practical importance of the 

problem (addition energies control the threshold poten- 
tials for one-electron transistors), experiments done in 
the low-density regime have revealed intriguing patterns 
[lit I which are suggestive of Wigner-like correlations and 



have never been fully explained. This kind of electronic 
structure calculations have been, so far, only accessible 
to wavefunction methods such as Quantum Monte Carlo 
[3] and configuration interaction (only small dots) [13] ■ 
The KS approach is, in this context, only useful in the 
moderately correlated regime [l8[ . Moreover, similarly to 
recently proposed first-order density-matrix energy func- 
tional [19|, the functional Ke^^[p] yields by definition 
the exact dissociation limit of the H2 molecule, a typical 
case in which restricted KS calculations fail. The more 
challenging case of moderately stretched H chains will be 
tested in future work. We believe that our approach will 
prove more effective than traditional KS for these kind 
of problems [j^ ]. contingent on the development of an 
efficient algorithm to routinely solve the SCE equations. 
A first step in this direction is to reformulate V^e^^[p] as 
a generalized mass transportation problem [20| . 

P.G.G. was supported by ANR (07-BLAN-0272), and 
G.V. was supported by DOE Grant No. DE-FG02- 
05ER46203. P.G.G. acknowledges discussions with K. 
Burke and A. Savin. 



[2: 

[4: 

[6: 
[7: 

[8 
[9 

[m 



[11 
[12: 

[is: 

[14 



[is: 
[16: 

[17 

[is: 

[19 



p. Hohenberg and W. Kohn, Phys. Rev. 136, B 864 
(1964). 

W. Kohn and L. J. Sham, Phys. Rev. A 140, 1133 (1965). 
M. Levy, Proc. Natl. Acad. Sci. U.S.A. 76, 6062 (1979). 
M. Seidl, Phys. Rev. A 60, 4387 (1999). 
M. Seidl, J. P. Perdew, and M. Levy, Phys. Rev. A 59, 
51 (1999). 

M. Seidl, J. P. Perdew, and S. Kurth, Phys. Rev. A 62, 
012502 (2000). 

M. Seidl, P. Gori-Giorgi, and A. Savin, Phys. Rev. A 75, 
042511 (2007). 

P. Gori-Giorgi, M. Seidl, and A. Savin, Phys. Chem. 
Chem. Phys. 10, 3440 (2008). 

P. Gori-Giorgi, G. Vignale, and M. Seidl, J. Chem. The- 
ory Comput. 5, 743 (2009). 

G. F. Giuliani and G. Vignale, Quantum Theory of the 
Electron Liquid (Cambridge University Press, New York, 
2005). 

M. Taut, Phys. A: Math. Gen. 27, 1045 (1994). 

C. Attaccahte, S. Moroni, P. Gori-Giorgi, and G. B. 

Bachelet, Phys. Rev. Lett. 88, 256601 (2002). 

R. C. Ashoori, Nature 379, 413 (1996). 

S. Tarucha, D. G. Austing, T. Honda, R. J. van der 

Hage, and L. P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 

(1996). 

N. B. Zhitenev, R. C. Ashoori, L. N. Pfeiffer, and K. W. 

West, Phys. Rev. Lett. 79, 2308 (1997). 

A. Ghosal, A. D. Guclu, C. J. Umrigar, D. UUmo, and 

H. U. Baranger, Nature Phys. 2, 336 (2006). 

M. Rontani, C. Cavazzoni, D. Bellucci, and G. Goldoni, 
J. Chem. Phys. 124, 124102 (2006). 

S. M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 
1283 (2002). 

O. Gritsenko, K. Pernal, and E. J. Baerends, J. Chem. 
Phys. 122, 204102 (2005). 



5 



[20] G. Buttazzo, L. De Pascale, and P. Gori-Giorgi, in prepa- 
ration. 

[21] The good performance of SCE-LDA also for large uj is due 
to the exactness of LDA for the non-interacting kinetic 
energy in a 2D harmonic trap, see M. Brack and B. van 



Zyl, Phys. Rev. Lett. 86, 1574 (2001) 
[22] Some similar formalism has been recently and in- 
dependently discussed by Z.-F. Liu and K. Burke, 
iarXiv:0907.2736t 



